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Abstract. A sample of 427 gamma-ray bursts (GRBs), measured by the 
RHESSI satellite, is studied statistically to determine the number of GRB 
groups. Previous studies based on the BATSE Catalog and recently on the 
Swift data claim the existence of an intermediate GRB group, besides the long 
and short groups. Using only the GRBs durations Tgo and x 2 or F-test, we have 
not found any statistically significant intermediate group. However, the maxi- 
mum likelihood ratio test, one-dimensional as well as two-dimensional hardness 
vs. Tgo plane, reveal the reality of an intermediate group. Hence, the existence 
of this group follows not only from the BATSE and Swift datasets, but also 
from the RHESSI results. 
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1. INTRODUCTION 

Originally it was found (BATSE data: Kouveliotou et al. 1993; Konus-Wind 
data: Aptekar et al. 1998) that two GRB classes exist: a short one of duration 
< 2 s and a long one > 2 s. They have different celestial distributions (Balazs et 
al. 1998, 1999; Meszaros et al. 2000a,b; Litvin at el. 2001; Meszaros & Stocek 
2003; Vavrek et al. 2004, 2008), and they contain two different central engines 
(Balazs et al. 2003, 2004; Fox et al. 2005). Some articles deal or even point to 
the existence of three GRB classes in the BATSE database (Horvath 1998, 2002; 
Mukherjee et al. 1998; Belousova et al. 1999; Balastegui et al. 2001; Horvath 
et al. 2004, 2006; Chattopadhyay et al. 2007). Some authors conclude that the 
third class (of intermediate duration), observed by BATSE, is a bias caused by an 
instrumental effect, placing the intermediate group as a separate source population 
in doubt (Hakkila et al. 2000, 2004; Rajaniemi et al. 2002). On the other hand, 
recent papers by Horvath et al. (2008), Horvath (2009), and Huja et al. (2009) 
indicate that there is a statistically significant intermediate group in the Swift and 
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BeppoSAX datasets. 

The purpose of this paper is to investigate the number of GRB groups in a 
dataset provided by the RHESSI satellite^] which consists of 427 GRBs covering 
the period of 2002-2008. This issue is examined in more detail in the paper Ripa 
et al. (2009). We analyzed both one-dimensional distribution of GRB durations 
and two-dimensional plane of hardness ratio versus duration. To determine the 
number of GRB groups, standard statistical tests described by Trumpler & Weaver 
(1953), Press et al. (1992) and Zey et al. (2006) were used. 

2. DATA SAMPLE 

The Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is a NASA 
Small Explorer satellite designed to study hard X-rays and gamma-rays from solar 
flares (Lin et al. 2002). Its spectrometer consists of nine Ge detectors (Smith et 
al. 2002) which are only lightly shielded, ensuring that RHESSI is also useful for 
detecting non-solar photons from any direction (Smith et al. 2003). The energy 
range for GRB detection extends from about 30 keV to 17 MeV. The effective area 
is around 150 cm 2 (Wigger et al. 2006). We used the RHESSI GRB Lisi0 and the 
Cosmic Burst Lis£| to find 487 GRBs between 2002 February 14 and 2008 April 
25. There is no automatic search routine, only if there is a message from any other 
instrument of the Interplanetary Network 0, the RHESSI data are searched for a 
GRB signal. 

We chose a subset of 427 events with signal/noise ratio higher than 6 and 
derived count flux curves (with a time resolution higher than 10 % of the burst 
duration for the vast majority of our entire dataset) and count fluences from the 
rear detector segments (except number R2) of the spectrometer (Smith et al. 2002) 
in the energy band from 25 keV to 1.5 MeV. 

3. DURATION DISTRIBUTION 

First, we studied one-dimensional duration distribution of Tgo (Kouvcliotou et 
al. 1993; Fishman et al. 1994). The histogram consists of 19 bins on a decimal 
logarithmic scale (Figures 1 and 2). The uncertainty of Tgo consists of two com- 
ponents. One is given by the count fluence uncertainty during Tgo (5t s ), which is 
given by Poissonian noise, and the second one is the time resolution of derived light 
curves (St res ). The total Tgo uncertainty 5t was calculated to be 5t = \/St^ + St^ cs . 
We follow the method used by Horvath (1998) and fitted one, two and three log- 
normal functions and used the \ 2 test to evaluate these fits. The number of GRBs 
per bin is at least 5 (except the last bin). The fit with one log-normal function is 
highly unacceptable, because \ 2 — 157 for 17 degrees of freedom (dof). The fit with 
two log- normal functions (dof — 14) gives \ 2 — 19.13 which imply a goodness-of-fit 
(gof) of 12%. The fit with three log-normal functions (dof — 11) gives x 2 — 10.30 
which imply gof = 41%. The assumption of two groups being represented by two 
log-normal fits is acceptable, the fit with three log-normal functions even more. 
The question is whether the improvement in \ 2 is statistically significant. To an- 
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Fig. 1. Duration distribution of the 427 RHESSI GRBs (the number of bins is 19). 
The best fit with two log-normal functions (left panel) and with three log-normal func- 
tions (right panel). The error bars are standard deviations of the number of GRBs per 
bin for ten different simulated duration distributions described in the text. 

swer this question, we used the F-test, as described by Band et al. (1997). The 
F-test gives a probability of 6.9% of the improvement in \ 2 being accidental. This 
value is remarkably low, but not enough to reject the hypothesis that the two 
log-normal functions are sufficient to describe the observed duration distribution. 
To determine how the Tgo uncertainties affect our result, we randomly selected 
one half of the bursts and shifted their durations by the full amount of their un- 
certainties to lower values and the second half to higher values. We then compiled 
a histogram and recalculated the best-fit model parameters, \ 2 an d F-test. Only 
three cases out of ten performed simulations gave F-test probability lower than 
5% (namely 0.5, 2.0 and 3.5%). Therefore, on average, the improvement in \ 2 is 
insignificant and we cannot proclaim the acceptance of the three groups of GRBs 
by the \ 2 method. 

Since the number of GRBs is just slightly higher than 5 for many bins, we 
also used the maximum likelihood method (see in Horvath 2002) to fit two and 
three log-normal functions. The difference of the logarithms of the likelihoods 
AlnL = 9.2 should be distributed as a half of the x 2 f°r 3 degrees of freedom. 
From this value we infer that the introduction of a third group is statistically 
significant on the 0.036% level. To estimate, how the Tgo uncertainties affect our 
result, we again generated ten different datasets randomly changing the durations 
by the full amount of their uncertainties. All ten simulations give probabilities 
which are in favour that introducing of the third group is accidental, much lower 
than 5% (in fact lower than 1%). Thus, the hypothesis of introducing a third 
group is acceptable when using a maximum likelihood fit. 

4. HARDNESS RATIO VERSUS DURATION 

A two-dimensional scatter plot (duration Tgo vs. hardness ratio) is shown in 
Figure 2. The hardness ratio is defined as the ratio of two count fluences F in two 
different energy bands integrated over the time interval Tgo- We used the energy 
bands (25-120) keV and (120-1500) keV, i.e., H = ^120-1500/^25-120- Employing 
the maximum likelihood method (see Horvath et al. 2004, 2006 and references 
therein), we fit two and three bivariate log- normal functions. The difference in 
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Fig. 2. Hardness ratio vs. Tgo of the RHESSI GRBs with the best fit of two and three 
bivariate log-normal functions. 

the logarithms of the likelihoods Alni = 10.9 should be distributed as a half of 
the x 2 for 6 degrees of freedom (Horvath et al. 2006) From this value we infer 
that the introduction of a third group is statistically significant on the 0.13% level. 
To estimate, how the uncertainties of GRB durations and hardness ratios effect 
our result, we generated ten different datasets randomly changed in durations 
and hardness ratios by the full amount of their uncertainties and recalculated the 
likelihoods. All ten simulations infer probabilities, of introducing the third group 
being accidental, lower than 5%. Thus, the hypothesis of introducing the third 
group is again acceptable. 

5. DISCUSSION AND CONCLUSIONS 

The existence of an intermediate class from the RHESSI Tgo distribution in 
not confirmed using the \ 2 method. However, the maximum likelihood ratio test 
on the same data reveals that the introduction of the third class is statistically 
significant. The hardness ratio versus duration plot for the RHESSI sample also 
demonstrates the existence of the third class. The typical durations are similar to 
those obtained with the one-dimensional analysis. The typical durations found for 
BATSE (Horvath et al. 2006) are roughly a factor of 2 longer than for RHESSI, but 
consistent for all three classes. The shorter durations of the RHESSI GRBs com- 
pared to the BATSE GRBs can be understood in the following way: for RHESSI, 
which is practically unshielded, the background is high (minimum around 1000 
counts per second in the (25-1500) keV band) and varies by up to a factor of 
3. Additionally, the RHESSI sensitivity declines rapidly below w 50 keV. Weak 
GRBs (in the sense of counts per second) and soft GRBs are not so well observed 
by RHESSI. Since GRBs tend to be softer and weaker with time, they rapidly fall 
bellow the RHESSI detection limit, resulting in a shorter duration being inferred. 

In the RHESSI dataset we confirm the BATSE result that short GRBs are 
on average harder than the long GRBs. The hardness of the intermediate class 
found for the RHESSI data is similar to that of the short GRBs. This is surprising 
since the intermediate class in the BATSE data was found to be the softest. This 
discrepancy might by explained by the different definitions of the hardness. The 
hardness H for the RHESSI data is defined as H = -F120-1500/-F25-1207 whereas for 
the BATSE data H = -F100-320/-F50-100, where the numbers denote energy in keV. 
This means that hardnesses measure different behavior of bursts. The situation 
differs even more significantly if we compare hardnesses in the Swift and RHESSI 
databases, because the Swift hardnesses are defined as H = -F100-150/-F50-100 and 
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H = ^50-100/^25-50 (Horvath et al. 2008; Sakamoto et al. 2008). 
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